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ABSTRACT 



The development of the three dimensional, laminar velocity profile 
in the entrance length of a rectangular duct is investigated. The solution 
to this hydrodynamic problem is obtained from the full, incompressible 
Navier-Stokes equations and the continuity equation, in finite difference 
form, on the digital computer employing the computational method of 
Chorin (l) . The solution yields the hydrodynamic velocities U, V, and 



W and the friction factor as a function of the distance from the entrance. 
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TABLE OF SYMBOLS 



English Letter 



A 


Cross sectional area of the duct 


d 


4 A 

Hydraulic diameter , 


D 


Perimeter of the duct 


F 


Friction factor 


P 


Pressure 


P 


Dimensionless pressure 


t 


Time 


T 


Dimensionless time 


u o 


Axial reference velocity 


u 


, Axial velocity 


u 


Mean axial velocity 


U 


Dimensionless axial velocity, u/u Q 


U 


Mean dimensionless axial velocity, u/u Q = 1 


V 


Crosswise velocity in y direction 


V 


Dimensionless crosswise velocity in y direction 


w 


Crosswise velocity in z direction 


W 


Dimensionless crosswise velocity in z direction 


X 


Axial normal coordinate 


X 


Dimensionless length, -|j- 


X 

e 


Dimensionless entrance region length 


y 


Crosswise normal coordinate 


y w 


Value of y at the sidewall boundary 
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Y Dimensionless normal coordinate 

z Crosswise normal coordinate 

z Value of z at the sidewall boundary 

w 

Z Dimensionless normal coordinate 



Greek Letter 



o Artificial compressibility 

A Crosswise coordinate, y or z 

^ Fluid density, artificial density 

T Shear stress 

yU Fluid viscosity 

V Kinematic viscosity 



Nondimensional Groups 

/°^ Reynolds number , R e 
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I. INTRODUCTION 



The design of compact heat exchangers for use in gas turbine 
regenerators , among other applications , entails the maximum use of 
the high heat transfer rates attainable in the entrance region. In this 
region the fluid is undergoing rapid transition from its uniform profile at 
the entrance to its fully established value farther downstream. A knowl- 
edge of the temperature profiles in the region is required to obtain the heat 
transfer rates. The temperature profile develops simultaneously with the 
velocity profile and is dependent upon it. The Graetz solution assumes a 
fully established velocity profile and thus eliminates most of the complexity 
of the hydrodynamic problem. This solution, however, fails to yield ac- 
curate results in problems involving large entrance region to overall length 
ratios . 

Various approximate methods have been devised for two dimensional 
channels, such as a circular pipe or parallel plate channels, to determine 
the development of velocity profiles. Until recently these solutions have 
invariably made boundary layer assumptions. Schlichting (5) obtained a 
solution for the flow between a pair of infinite parallel flat plates. He 
used two asymptotic series solutions, one based on Blasius' solution of 
the boundary layer development expanded in the downstream direction , 
and the other based on the Hagen-Poiseuille solution of a parabolic 
velocity distribution, expanded in the upstream direction. He then 
joined these two solutions. 
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Recently Wang and Longwell (6) have employed the full momentum 
equations and the introduction of a stream function to obtain a solution 
for the velocity development between parallel flat plates 0 The resulting 
equations were cast in finite difference form and the solution effected 
numerically. To date, no exact, three dimensional solutions for develop- 
ing velocity profiles have been obtained. 

The purpose of this analysis is to present an exact, finite- 
difference solution to the full incompressible Navier-Stokes equations 
for the hydrodynamic entrance region of a square duct. This solution is 
not nearly as simple as approximate methods. Its value, however, lies 
in its potential accuracy and the fact that it is an exact solution on 
which the credibility of various approximate solutions may be based. 
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II. ANALYSIS 



A„ Governing Equations 

The governing equations are the momentum (Navier-Stokes) and the 
continuity equations. The complexity of the equations requires the fol- 
lowing assumptions concerning the flow: 

1 . steady and three dimensional 
2 o laminar 
3. incompressible 

Through these assumptions, the energy and momentum equations become 
uncoupled, and may be expressed as follows: 

Case I: 



Momentum: 















(i) 












■ u 



Lw * 



o- 



^ 3+ d 



cf -1 



( 2 ) 



Auj + n <W . = - -L AB. \ -\>\ + 4-4^- 



(3) 



Continuity: 

^ *A. 4 sL-i£ 4 _ 



(4) 
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B. 



Boundary Conditions 



The boundary conditions associated with the equations for a square 
duct with a uniform velocity at the entrance are: 

1 . At the entrance ( x = 0) 
v = w = 0 



U = U Q 

2 0 At the walls (x = 0, y = y w , z = 0, z=z w ); 

u = v = w = 0 
Non-dimensionalization 



Equations (l) , (2), (3) and (4) may be written in dimensionless form 
through the introduction of the following dimensionless variables: 



u-tt. v=^ 



v _ J<_ Y--M_ 7- = jL T- 

a 1 a c a 

The non-dimensional equations are then: 



(5) 



Momentum: 



^ -t- R 



V ^Y w ^2 



^X ^x 2 ^Y 1 ^i 3 - 



(6) 



IY R 4 u w +v ^ w i 



^Y ^X 2 ^Y 2 ^z 2 



(7) 



^ w 





^2J ^2 ^Y 2- ^2 



( 8 ) 



Continuity: 

4X ^Y ^2 " 
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(9) 



D . Friction Factor 



To obtain an equation for the friction factor, the relation between 
the shear stress and the flow variables must first be introduced: 



a- u 



or r. 

y-"° 7 ^ 



y 



= o 



do) 



Since the duct is symmetrical, either derivative, or * may be 
used in equation (10). By introducing a new variable Y|, where T^= y or ^ , 
equation (10) may be written: 



r= 






( 11 ) 



The friction factor, F, may be related to the shear stress by the equation. 



f - JLL 

r f ui x 



(12) 



By substituting equation (11) into (12) the desired equation for the friction 



factor is obtained: 



F = ^ 



4 






f U 1 dV\ 






(13) 



By defining a new dimensionless variable, 7^- , equation (13) may 

be written in dimensionless form 



F s 2~ 

±y\ 



(14) 



y[- o 
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III. METHOD OF SOLUTION 



A . Technique 

The recently formulated technique of Chorin (l) was adopted for the 
solution of the governing equations . The principle of the method of solution 
suggested by Chorin lies in the introduction of an artificial compressibility, 
£ c into the equations of motion in such a way that the final results do not 
depend on £ . The solution is in essence a relaxation method. In employ- 
ing this method, it is necessary to introduce a modified set of governing 
equations: 



Momentum: 






R. 







(15) 












lY ax 1 






( 






f ^ £W 

^Y 5 U 1, 



Continuity: 

^ ^ 4. 4. ^ v_ j. ^ w - q 

^Y 4 % 



State: 

p= e A 



(16) 



(17) 



(18) 



(19) 
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In this set of auxiliary equations (f is an artificial density, $ 
is an artificial compressibility, and P= is an artificial equation 
of state. 

Equations (15), (16), (17) and (18) may be put into finite difference 
form by introducing the following finite difference approximations: 

uQax,y, 2 ,t) _ u(x t )-u(x') 

AX " “ AX 



U(x+-aX j Y ) z,t)i-U(x.-AX,Y,2- 1 t) - U(x,Y,2: J+ at) -U^.Y.ijT-AT) 
AX 2 ‘ AX a 



- + U(X~)-U(T + )-U(T~) 

AX 5 

The set of modified equations put into finite difference form and 
solved for the appropriate variables yields the following computing 
equations: 



U (T> 



± 



1 + 2AT+ 3 AT 3AT 

ax 1 ay 1 az 1 



-R e AT 
AX 



uV) - uV) 



- ^AlTii(r)v(y + )~u(y')v(Y')]- ^l[u(z T )w^)-^w(2 

AY L J A2 . 



+ DaT 
AX 1 



a. DAT 
A^ 



u(x + )+ulx')-u(T*) 



DAT 

AY Z 



u(Y>u(v-)-u(r‘) 






u(2 + )t u(2') - ur-)J-f^-[e(X + )-C(X-)] +■ 0 (T-) 



(20) 
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V(T> 



1 



r 



, . 2AT 3 AT + 2AI i 
1+ AX’ + AY 3 A ? 1 



juxO v(\ r ) - u(x-) V(x') 



- mi 

AY 



\jX V) - V l (Y)1 - %£ vJ(2 + )V(z + ) - W(r) v(z) 

bt 



a>T 

ax z L 

sat 
LI'l 



V(x^v(x')-V(r-) 

v(H + )+v(r)-v(T‘) 






-il 

SAY 



v(Y + )+v(r)-vi 



f(Y*) 




~\ 



J 



( 21 ) 



, + iAT ^AT^aAl 
^“AX 1 "ST‘ ai^ 



. R £ AT 



MI 

AX 



U(x + )W(x + )-^( x ')^CX‘) 



AY 

.aAT 

AX 3 

•^Al 

A 2 1 




v(y + )vj(y + )- 
W(x + )tW(x*)-W(Tj 
w(z + )+W(z ) ~w(t ) 



MI 

A* 



wV)^) 1 



:at 
AY 3 L' 



wM+w(V-)-w(r) 



j 

"\ 



AT, 

4 A 2 . 



f(z f )-?fe") +wfr) 



(22) 



^ ( T t) = -g j^u(x+) - U (X-)]- [ v ( yt ) ■ ' W Y ‘) - f^[w(z + )- wfcjJ^CM 



( 23 ) 
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B . Boundary Conditions 

In dimensionless form the boundary conditions for equations (15), 
(16) , (17) ; and (18) are: 

1 . At the entrance: (X = 0): 

V = W = 0 
U = 1 

2 . At the walls (Y = 0 , 1 Z = 0 , 1). 

W = V = U = 0 

3 o At X = X entrance length: 

W = V = 0 



C . Stability 

With the introduction of the modified equations (15) , (16) , (17) , 
(18) , and (19) an artificial Mach number is also introduced. 



Chorin (l) has stated that, in order to insure stability, the flow 
Mach number must be kept less than one. An additional requirement that 
must be met for stability of the set is that: 



^ ^ fully established 





(field) 
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D. 



Calculation of F 



In order to calculate the friction factor, it is necessary to 
introduce a forward difference approximation as follows: 

uCn+AYQ-U(y^ = u(n + )-U(^ 






Using this approximation, equation (14) may be written: 



8_ UCV) 

^ Re UY[ 



■q = o 



(24) 



Since equation (14) is applied at = 0, U (fj ) is identically zero and 

is eliminated from the finite difference expression. 

Equation (24) may be used to obtain an average value of the friction 

AU I 

factor. An average value of | - _ Q may *- >e determined at each X-wise 
location by summing U (f( + ) at each grid point on the boundary and 
dividing by the number of points used in the summation. The average 
value of the slope of the velocity profile may now be used to calculate 
an average friction factor as follows: 

6 



F = 



^ W+ 1 



(25) 



W AY\ F U 

= number of points used 

Uwti = velocity at one grid point from the side wall boundary 



The grid system and summation points are shown in Figure 1 . 



E . Numerical Procedure 

A finite difference solution was obtained, using the IBM system 360 
digital computer. Equations (20), (21), (22), and (23) involve approximating 
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a solution for the entire flow field at time equal to T and T and then 
calculating a new approximation for time equal to T + The computed 
values of U , V, and W are tested for continuity in equation (23) , and 
if C(T + ) *t(T‘) does not equal zero, then a steady state solution 
has not yet been achieved. T is then set equal to T + and T~ equal to T. 
Then the process is repeated until a steady solution is reached. That 
is until = O o This sequence serves to relax the set of equations. 

Solution also requires that be approximated at the grid points 

on the side wall boundary since equation (23) cannot be used at these 
points. To apply equation (23) , one would have to impose the condition 
that the flow parameters are identical on either side of the channel walls 
This assumption, however, is unrealistic. Instead, the following approxi- 
mation is introduced at the side wall boundary: 

-ME V( V) + f(r-) 

Using the stability requirements, the solution to equations (20), (21), 

(22), and (23) was obtained using the following values of grid parameters: 

AX=.00lM AT" - .oooooi AY*A2=.I ga.lfc 
The value of At was obtained using the entrance length of .0328, 
computed by McComas (4) , for a hydraulic diameter of 1.0 and a 
Reynolds number of 1.0. The fully established velocity profile at this 
entrance length was taken from Knudsen (2) . In this series solution for 
the fully established value of U, the pressure gradient was replaced 

by C k) , where F is the fully established factor, also taken from 

Knudsen (2) . 
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The value of F computed by Knudsen (56.24) agrees well with the 
value suggested by Lundgren (3) (56.908) and MoComas (4) (56.908). 

The computer program written in Fortran IV is included in the 
appendix. 
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III. RESULTS 



Various combinations of AX, AY, Al, AT, and £ were used in the 

solution, all of which satisfied the stability requirements set forth by 

Chorin (l) . It was also found, however, that unless the ratio 

was of the same order of magnitude as the other ratios of this type, the 

solution would diverge. The values of AX, A YjAi, AT, and £ given 

previously were found to be the most rapidly convergent values tried. 

Profiles were obtained using these values and Re equal to 1„0. They 

are plotted in Figure 2. The development, of the centerline velocity is 

shown in Figure 3 and Table I. Contrary to expectations, it was found 

that for any value of entrance length , the flow would not become fully 

established before it reached the X imposed upon it. This was true 

e 

even though the X g was much larger than the value computed by McComas 
(.0328) . In order to test the validity of the method, a solution for fully 
established flow in the duct was obtained using the fully established 
profile as boundary conditions at X = 0 and X = X entrance length. In 
this case, the profile at every X-wise location should be the same. 

Since the X-wise station X= Xe is the last station to be effected 
by the relaxation technique, the profile at X = Xe/ a. is plotted. For 200 
iterations the solution has not yet converged, and the probable error based 
on a calculation of Vfae ) on the centerline is 15 4 per cent. After 
500 iterations, the probable error has decreased to 0.7 per cent, and after 
800 iterations, the probable error is reduced to 0.1 per cent. One can 
see that as the number of iterations is increased, the accuracy of the 
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solution is increased until a very high degree of accuracy (99.93 per cent) 
is obtained at 1500 iterations. This is shown in Table II, Figure 4 and 
Figure 5. 

The average friction factor obtained as a function of X in the entrance 
length is shown in Figure 6 and Table III. The first four to six computed 
points were in error due to the finite value of aY and A 2- , which yield 

a finite value of shear stress at the entrance. In actuality the shear stress 
is infinite at the entrance, due to the discontinuity in the velocity at 
Y or 1 ~ O • In the limit as approaches zero, equation (25) does, in 

fact, approach infinity. As the entrance length is approached, the value 
of F asymptotically approaches a value of 54. This value is in very good 
agreement with the fully established value that was introduced into the 
boundary conditions . It also is in good agreement with the value suggested 
by McComas (4) and Lundgren (3) . 
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V. CONCLUSIONS 



As a result of this analysis , it may be concluded that: 

A solution using a finite difference approximation, and relaxation 
type of scheme suggested by Chorin (1) , leads to a convergent and 
accurate solution albeit with substantially more computation than is 
inherent to an explicit scheme. The accuracy of the solution depends on 
the number of iterations made, the size of the finite-difference variables 
AX i , A and hT , and the value of the entrance length imposed 
upon the flow. This dependence of the solution upon the length of the 
entrance region is not what was expected and should be the subject of 
further investigation. The solution obtained in this analysis was based 
on a Reynolds number of 1.0. As one can see, the computing equations 
contain Reynolds number as a parameter and therefore depend upon it. 
Before a solution of this type may be used in the calculation of heat 
transfer in an entrance region, it would be necessary to investigate the 
behavior of this method at Reynolds numbers that are more typical of 
actual flow problems (R e = 10 ). A further limitation to the solution 
obtained in this analysis is that of computer time and storage space. 

This might be overcome by making use of the symmetry of the duct with 
a suitable modification of the boundary conditions. 
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APPENDIX A 

SAMPLE FORTRAN IV PROGRAM 
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C THIS PftWRAM EVALUATES THE FULL TNCOMPP FS S I BL E NAVTFP-STOKES 

L AMD CONTINUITY EQUATIONS IN THE ENTRANCE R pOlON OF A SQIIARF DUCT 

c VIFLOINi; THF HYDE ODYNAM I C VAP TABLE S U,V,W AND P AT FVFRv GRID 

C POINT »N THF ENTRANCE LFNOTH. 

C 

c 

c 

c 

DIMENSION H<?l,ll.llt3I.V<21«U,ll,3).W(21til.ll,3l,P(?!,ll,ll,3»' 
READ! 5, 1 ) I I , JJ ,KK, MM, NN 
1 FORMAT ( S I 6 I 

RFAD< 5, 2 IPE. A.DT ,OX,DY,DZ ,OFl TA 
? FORMAT ( 7E1L-. 6) 

N = 1 

P l = 3. lA15F?65358c:7 

F=?4./RF*l./(l.-19?./PI**5*(TANH<Pl/2. I - 1 « /3. **5 * T ANH ( 3 . * P I /? . MI ♦ 
I I 1 = 11-1 

DO 3 i = 1 » I I 
DO 3 J= 1 . MM 
DP 3 K = 1 ♦ N N 

u ( 1 1 j . k , 1 1 =o . n 

U I I.J.kfpl— I • u 1 
U< I ,J,K ,3) =.' ■ .( 

V ( I * J tK » 1 ) =(, #L 
V( I*J»K»2)=C.O 

V ( I , J . K , 3 I =f . C. 

W ( I » J » K « 1 I =0.0 
W ( I » J ♦ K ♦ 2 ) =f *1. 
w( I . J . K , 3 > => > • C 
P f I * J * k « 1 I = * C. 

Pd.J.K . 2 I =« .( 

I * J * K . 3 ) =t .{ 

3 CONTINUE 
DC 11 M= 1 , 3 
on u j=2.jj 
no 11 K=2fKK 

Y=-. 5+U-l I ✓ 1 C . 

Z=-. 5-MK-l I /1( . 

U(1.J.K,M)=1.0 

U( II . J,K , Ml = ?. /P I **l*F*PF*( f 1 . - COSHL P l*Y I /COSHl PI/?. ) l*CDS( P?» 7 ) 
*-l./27.*< l.-CnSH(PI*3.*Y)/COSH(PI*3./?. U*COS( Pl*3.«7 ) 
v ♦ 1./ 12 5.* < 1. -t.OSH< P 1*5. *Y) /COSH ( P[*5./?.n*COS<PT<'*.*-7»» 

11 CONTINUE 

10 CONTINUE 

DO 5 1=2.111 

DO 5 J= 1 . MM 
00 5 K = 1 « NN 

AA=1. ♦2.*DT/DX**2+2. *DT/DY**2+2. +DT/D7* *? 

IFIU.EO. 1 IPO TO 13 
IP(k.E0.1 IPO TO 1? 

IFIJ. tQ.il IPO TO 13 

IFIK. FO.ll IPO TO 13 
C 

C 

U( I fJ.K f 3l=(-RF^nT/DX*(U( I«-l.J.Kfll*<«2-U( I-l.J.K.l I’M?) 

*-R E*DT /DV* (IJII , J* 1 , K , l I*V( I.J + l , K, 1 I -U( I,J-1,K,1I*V( f.J-l.K.l»> 

*+2. +DT/0X**2*( 0( I«-l. J.K» 1 )+U( I — 1 . JtK.l *-ur I ,J,K, ?)) 

^♦2.*m/DYx *2*< U< I,J*1,K.1I+IIII,J-1,K,1HI(!,J,K,2M 

*-PF*DT/D7» (IK I , J.K 4 1 , 1 I *W(I , J.K + l , 1 I -III I f J ,K-1 , 1 I ^W( I , J.K-l , t * » 

♦ ♦2.*OT/D7’i t *2*(U( I .J.K+l.l I •*- IJ ( I * J ♦ K - 1 .1 >-ll( I ♦ J » K , 2 ) I 
»-DT/DX*l. /DEL TA*( P ( I ♦ 1 . J.K , 11 -P( 1-1 , J,k , 1 | I 4l»( f , J, K, 2 II* 1. / AA 

C 

c 

V(I,J,K,3I = ( - P F ♦ D T / D Y * ( V( l.J+l.K.l I v + ?-V( I ♦ .1- 1 * K , ] 1**2) 

X‘-RF*DT/nx* I U( I 4 i , J f K , 1 ) *VU + 1 . J »K , 1 I -ll( 1-1 , J.k , l ) *V( T-l . J ,K , 1 » | 

»-RF*DT/D/M W( J ♦.) ,K 4 1 , 1 I *V( T , J ,K 4 1 , 1 » - w ( I , .1 ,K-1 , M>iV( I . .1. k - J , 1 I) 

4 2 . *0 T / D X » * 2 * ( Y (I 4 1 , J , K , 1 I 4 v ( I - 1 . J . K ♦ 1 I - \/ < I . J . K , ? II 

♦ 42, =D1 / DY442M VI I , J4l ,k , 1 ) 4 V( I «.l-l ,k , l I -V( I , J,k , 71 I 

♦ 4?.v{ ( T/n/4 = 24(V( I , ,k 4l , l I 4 v ( I ..I.K- ! . 1 1 — V ( I » J • k » 2 I » 
*-DT/0YM./DFLTA4<r>(I f j4l,K,ll-P<t,J-l,K.l*l4V(I,J,K,?)I*l./AA 

(. 

r. 

W( L, J.K, 31 =(-RE«DT / D?*( W( ItJ.K+1 * 1 )**?-<*!( T «J»K — 
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* -RF*l>T/0**(U< I 4 1 , J.K , 1 )*WI I + 1, J.K, M -l 1(1-1 , J,* *1 )*Wf I-l, J , * . 1 » I 

*-«t*OT/D¥ti<V{ | ,-j+l , 1 )*W{ l , j-M , * , f ) -\t( T , K tl I * w( 1,1-’ . k , i || 

M?.^m/DXvv;< I M I ♦ 1 , J , K, l.» *-W( I - 1 ,J , K , 1 ) -W f I , J , K , ? ) > 

* *2. *01 /0vvw<2“( M| t ,)U ( K,lltW(l t J-l.K,ll-W(l,J,K,?ll 
•♦2.*CT/n7**2*(Kf I ,J.K*1,1 »4W( I » J , K - 1 ,1 )-Wf I , J, K . 2 ) ) 

«‘-DT/P2*l./r)ELTA*<P(It.J,Ki-l,l>-P<T,J,K-l.in«-W<i,J,K,2M*l./AA 

C 

C 

o I I , J , K , 3 ) ^-DT/OX* (IJ( I ♦ 1 , J ,K , 1 )-Ul 1-1 , J, «, 1 H 
*-DT /0Y*< v< I , J* 1 ,< . 1 l-v ( I . J-l , K , 1 M 
»-l)T/OZ*( W( I , J , K* 1 , 1 ) - W( I ♦ J ,K-1 , 1 ) ) 

**PU,J,K,?> 

1’ CONTINUE 

c 

IF I K. FO. I 1. AND. J.l T. 1 1. ANO. J.GT. 1) 

*P< I. J.li,3l = -3.*DT/P7*<W! I , J.1C, 1)) 

* + P< I , J, 11 . 21 

IF Ik. EG. 1, AND. J.OT. 1. AND. J.LT. II ♦ 

*P(I.J.l . 3 )=-"*. *DT /U7* I W( I ,J,2 ,1)> 

^ ♦ P f t.J.l , 2 ) 

IF ( J.EO.l I. AND. K. IT. Il.AND.K.GT. 1 ) 

*P( I , 11 »K, 3 ) =-3. *DT /l)Y*I VI I , 10 ,K , 1 M 
♦♦PI I » I 1 » K , 2 ) 

IFIJ.EO.l . AND. <. GT. 1. AND. K.L T. 1 l ) 

*P(I»l ♦K,3Is-3.*DT/nv*(viI,2 » K , 1 ) I 
♦♦Plltl ,*.2> 

5 CONTINUE 

DO 6 1=2.111 
on 6 J = 2 ♦ J J 

DO 6 K-2.KK 

UI I . J .K , 1 ) =Uf I .J ,K . 3 ) 

VI I , J ,K . 1 ) =V( I , J ,K ,3) 

* ( I »J»K»1)=W( I * J * K . 3 ) ' 

P< I . J.K. 1) =P(I . J.K ,3) 

U I J , J.K , 2) =01 I ♦ J » K. » 1 ) 

VI I , J.K . 2 » = V I I ,J ,K ,1 ) 

Wl ! , J ,K ,21 = W< I , J.K , 1 ) 

PI I «J.K,2I=P( I .J.K ,1 ) 

6 CONTINUE 
URITEIF.15IN 

1 FORMAT! 15 1 

IPIN.F0.20’ » GC! TO 15 

IFIN.EQ.5.,0 ) GO TO J6 
IFIN.E0.8v,; »C( TO 16 
IFIN.E0.11 C 1 Of' TO 16 
IF (N.PO. 1500 IGP TO 16 
N = N+ 1 

on to in 

16 CONTINUE 
1 A CONTINUE 
WR I TF I 6 , 7 » 

7 FORMA TI//T3.' POSITION', T27,»U»,T45,»V',T6 3. IW'.TRl.'P*, Too, .T‘,//) 

m = 3 

DO <t 1 = 1.1 I 

no 0 J= 1 , MM 

CO 9 K = 1 , NN 

PP = P( I , J.K .M) /DELTA 

TT6I 6 ,8 M.J,K,u(I.J.K,M»,V!I,J,K,M*,wit,J,K,M»,PP.N 

8 FORMAT (315 .4F1P.8, I5» 

9 CONTINUE 

138 FORMAT! ///T3. 'U-CFNTERL INE/U-MAX* ,//) 

WRJTF (6,138) 

DO 137 1=1,11 

UCUM = IJ( I .6,6,3 I /'M 21 , 6, 6.31 
WRITE (6. 139) Uf UM , I 

139 FORMAT I F 1 8 . 8 » 16) 

137 CONTINUE 

N = N4 1 

IFIM.LT. 153DIG0 TO 1C 

P.ETUR N 

END 
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GRID POINTS USED IN 
THE SUM FOR F 



FIGURE I 
GRID SYSTEM 
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FIGURE 2 

VELOCITY PROFILE DEVELOPMENT 




FIGURE 3 

DEVELOPMENT OF THE 
CENTERLINE VELOCITY 
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VELOCITY PROFILE DEVELOPMENT 




FOR FULLY ESTABLISHED FLOW 
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FRICTION FACTOR 



X 



TABLE I 

CENTERLINE VELOCITY 
DEVELOPMENT 




0 

.00328 
.00656 
.00984 
.01312 
. 0 I 6 4 0 
. 0 I 9 6 8 
.02296 
.02624 
.02952 
0 3 2 8 0 



.4815 
.5314 
.5816 
.6 3 21 
. 6 8 30 
. 7 3 43 
. 7 8 62 
0 3 86 
.8916 
.9 453 
1 . 0 
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TABLE II 

PROBABLE ERROR 
OF THE 
SOLUTION 



NUMBER OF 
ITERATIONS 



U 



% 



AT X*X, 



MAX 



(COMPUTED) 



U 



Vu 



AT X«X, 



MAX 



(ACTUAL) 



PROBABLE 
ERROR (%) 



2 


0 0 


. 8 4 


5 7 


5 


1 . 0 


15.4 




5 


0 0 


. 9 9 


2 8 


8 


1 . 0 


. 7 




8 


0 0 


. 9 9 


9 0 


8 


1 . 0 


. 0 


9 


1 1 


0 0 


. 9 9 


9 3 


2 


1 . 0 


. 0 


7 


1 5 


0 0 


. 9 9 


9 3 


3 


1 . 0 


. 0 


7 
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TABLE III 

CALCULATED FRICTION FACTOR 



X 



0 

.00164 
.0 032 8 
.0 049 2 
.0 065 6 
. 0 082 0 
. 0 098 4 
.01148 
. 0 13 12 
. 0 164 0 
.01968 
.02296 
.02624 
.02788 
.02952 
.03116 
. 0 328 0 



8 0 
7 8.1 
7 6.3 
7 4.3 
7 2.6 
7 0.17 
6 9.05 
6 7.3 
6 5.7 
• 4.16 
6 1.1 
5 8.7 
5 5.7 
5 4.82 
5 4.64 
5 4.4 
5 4.2 
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